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Abstract 

We develop a method for computing exact one-loop quantum corrections to the ener- 
gies of static classical backgrounds in renormalizable quantum field theories. We use 
a continuum density of states formalism to construct a regularized Casimir energy in 
terms of phase shifts and their Born approximations. This method unambiguously 
incorporates definite counterterms fixed in the standard way by physical renormal- 
ization conditions. The result is a robust computation that can be efficiently imple- 
mented both numerically and analytically. We carry out such calculations in models 
of bosons and fermions in one and three dimensions. 
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Chapter 1 
Introduction 



Given a static field configuration in a renormalizable quantum field theory, a natural 
question to ask is what its energy is. This question is especially relevant in the 
study of solitons, since they are local minima of the energy. Classically the energy is 
simple to calculate from the Lagrangian density of the theory. However, quantization 
of the theory introduces corrections to this classical result. Such corrections can be 
expanded as a power series in the coupling constants of the theory, which is equivalent 
to an expansion in h. 

In this work we will consider the leading quantum corrections to the classical en- 
ergies of field configurations in a variety of quantum field theories. Since the power of 
h counts the number of loops in the diagrammatic expansion of the energy, taking the 
leading correction will correspond to a one-loop calculation, in which we sum all one- 
loop diagrams with arbitrary numbers of insertions of the classical background field. 
This calculation is equivalent to summing the shifts in the \hui zero-point energies of 
all the small fluctuations modes in the presence of the background field. This sum 
is known as the Casimir energy. As usual in quantum field theory, the result of such 
a calculation diverges, and we must introduce divergent counterterms that are fixed 
through a finite set of renormalization conditions. The renormalization conditions 
define the theory in terms of physical quantities. Our challenge will be to implement 
such a calculation in a robust, efficient and unambiguous way. In particular, we must 
isolate the cancellation of divergent quantities without missing any finite contribu- 
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tions, and we must be sure that we have implemented the renormalization conditions 
faithfully. Merely "canceling the infinities" is clearly not sufficient when we want to 
regard a finite result as a physical prediction of a particular theory defined under 
fixed renormalization conditions. 

One application of this method is in coupling of the Higgs sector of the Standard 
Model to heavy quarks. If we imagine adjusting the Yukawa coupling for a quark dou- 
blet so that the quarks' mass becomes very large, we would expect that when they 
become sufficiently heavy, they effectively decouple from the theory. However, this 
cannot be the whole story, since if we simply removed a quark doublet from the Stan- 
dard Model, we would ruin anomaly cancellation. The resolution of this paradox Jl| 
requires that solitons in the Higgs field carry the quantum numbers of the decoupled 
fermions. This result suggests a picture in which heavy quarks are realized at small 
coupling as elementary fermion excitations, and at large coupling as Higgs solitons. 
In between these two limits, one might then expect to find a hybrid configuration, 
with the heavy quark tightly bound to a deformation in the Higgs field. To see if this 
picture is correct, we need to do a variational computation: we must find the field 
configuration of lowest energy that carries the heavy quark quantum numbers. This 
application highlights the importance of unambiguously fixing our renormalization 
conditions and avoiding finite errors in our energy computation. If our renormaliza- 
tion conditions were not fixed precisely, we would effectively be changing the theory as 
we moved from one background field to another, rendering the variational calculation 
meaningless. 

One could imagine building a Higgs configuration that brings a heavy quark level 
down from the mass of the quark closer to (or below) zero. The Higgs configuration 
itself would have an energy cost from the gradient and potential terms in the classical 
bosonic Hamiltonian. This cost could be balanced by the shift in the "valence" 
quark energy level, since the quark occupying this state moves to a lower energy. 
In particular, for strong quark Yukawa coupling (which means large quark mass), 
and small Higgs self-coupling (which means small Higgs mass), deforming the Higgs 
field would be favored. However, these two pieces alone do not form a self-consistent 
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semiclassical calculation. We have included the classical energy along with one part of 
the leading quantum correction to that energy, the shift in the valence level. We have 
no justification for ignoring the shifts in all the other quark levels, which contribute 
at the same order through their zero-point energies. Thus we are forced consider the 
full Casimir energy in this problem. If we do include the Casimir energy, we obtain 
a self-consistent semiclassical result. One may doubt the validity of the semiclassical 
approximation as the coupling gets large, of course, but we will obtain a result that 
is valid in a well-defined approximation. (In some models, it is also the exact result 
in a large- N limit ||.) 

Although we must compute the full Casimir energy, not just the valence contri- 
bution, it is still possible to construct self-consistent approximations to this quantity. 
Of course, these approximations further restrict the domain of validity of the com- 
putation. In particular, if the background field is slowly varying on the scale of the 
Compton wavelength of the quantum fluctuations, the derivative expansion becomes 
valid. Indeed, it has been used in models similar to the ones we will consider [|J. 
However, in both models of heavy quarks and other models we will want to explore, 
the scale at which the background field varies is precisely the Compton wavelength 
of the quantum fluctuations, so that all the terms in the derivative expansion will be 
about the same size, rendering the expansion unreliable. Our approach will be exact 
to one-loop order. 

Central to our technique will be to re-express the Casimir sum in the continuum 
in terms of phase shifts, using a formalism that originates with Schwinger's work on 
QED in the presence of strong fields [Q. We will consider only field configurations with 
some form of spherical symmetry, so that we can use a partial-wave decomposition. 
In each channel, the difference between the free and interacting density of states is 
related to the phase shift by 



which follows from imposing a boundary on the system and then sending the boundary 



p(k) = p (k) + 



1 dS(k) 




ir dk 
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to infinity, or by more formal S-matrix arguments. In any channel, for any k, the 
phase shift is a finite, well-defined, physical quantity. Its analytic structure is well 
understood in terms of Jost functions, and it can also be rigorously related to the 
Green's functions and S-matrix of the theory. The sum over continuum modes is then 
replaced by an integral over the density of states 



in which the bound states are still included explicitly. This mathematical artillery 
leads to several important properties of the phase shifts: 

• We can expand the phase shifts as a Born series in the strength of the poten- 
tial. This expansion is in exact correspondence with the expansion of the full 
propagator of the theory in terms of the free propagators connecting insertions 
of the potential. The Born expansion also has simple behavior at large k and 
becomes more and more accurate in this limit, which will enable us to use it as 
a regulator of ultraviolet divergences. 

• The phase shifts track level crossings, ensuring a consistent counting of modes 
as states become bound. This property allow us to avoid a serious problem we 
would encounter if we put the system in a box, because in that case we would 
have a hard time ensuring that we have kept the appropriate number of modes 
when computing the Casimir energy and the contribution of the counterterms. 
Missing even one mode in the Casimir sum will lead to a drastic change in the 
final answer; although such an error is small compared to the leading (divergent) 
behavior of the sum, the leading behavior is cancelled by the counterterms. The 
final answer is generically of the same order as a typical energy level. Our key 
tool here will be Levinson's theorem, which relates the phase shift at k — to 
the number of bound states. 

• The phase shifts and their Born approximations are simple, robust quantities 
amenable to both numerical and analytic calculation. We will see that they 




(1.2) 
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enable us to replace the cancellations of large quantities, the Casimir energy 
against the counterterms, by the much more manageable cancellation of exact 
phase shift against Born approximation. Many cases will in the end be tractable 
only numerically, but we will also find that numerical analysis will also shed light 
on analytic results by allowing us to continuously interpret between a trivial 
configuration and one that can be solved analytically. 

• The absence of boundaries will allow us to avoid spurious contributions from 
artificial boundary conditions. This property will be especially useful in models 
with fermions. 

• In principle, phase shifts can even be computed in fractional dimensions, since 
they obey a simple radial differential equation that depends analytically on 
the dimension of space. In particular, the first Born approximation the bosonic 
phase shifts can be analytically calculated in arbitrary dimensions terms of gen- 
eralized Bessel functions. This result agrees exactly with the result one would 
find for the corresponding tadpole graph evaluated in dimensional regulariza- 
tion. 

Much of this work originally appeared in Q. 
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Chapter 2 

<\y" theory in 1+1 dimensions 



We will begin with examples in 1 + 1 dimensions. The study of exactly soluble 1 + 1 
dimensional problems has yielded many insights into fundamental problems in field 
theory. Other 1+1 dimensional problems cannot be solved exactly, making it impor- 
tant to understand which properties of exact results will generalize to more generic 
cases and which are special to the exactly soluble cases. The renormalization process 
is simpler in one dimension than it is in three, since theories in lower dimensions are 
less divergent, but it still must be approached carefully. Indeed, we will find that the 
one dimensional problem actually contains additional subtleties not present in higher 
dimensions. 

2.1 Formalism 

We will consider a standard 4 theory with spontaneous symmetry breaking and a 
source J(x). The action is 

S[4>] = ^ / (^{d^f - ^(0 2 - l) 2 - J{x)4> + C^j d 2 x (2.1) 

where C ct is the counterterm Lagrangian. Our metric has g u = —g X x = 1- We 
have rescaled the field somewhat unconventionally in order to make explicit the 
correspondence between the powers of A and the powers of H, which we have set equal 



11 



to I. Classical terms will go as j, the one-loop terms we compute will go as A , and 
higher loops will contribute with higher powers of A. The mass of fluctuations around 
the trivial vacua (f>(x) = ±1 is m, and we define the potential U((f>) = ^((fi 2 — l) 2 . 

We consider a fixed field configuration (f>o(x). Initially we will assume that <f>o(x) = 
4> (—x), which restricts us to the topologically trivial sector of the theory. However, 
we will see that our method works equally well for the case of 4>q(x) = —<f>o(—x), so 
that in fact we can deal with configurations with any topology as long as U (0o) has 
reflection symmetry. We adjust the source so that </>o is a stationary point of the 
action, which means that J then solves the equation 

J(x) = ^-lm 2 ($-to). (2.2) 

We can always solve this equation for J, but there is no guarantee that the J we find 
will always correspond to a unique 0o, as we will see later. If O is a solution to the 
equations of motion, of course the source will be zero. 

We would like to consider the leading quantum correction to the classical energy 
of this configuration, which we can represent as the sum of the zero-point energies of 
the normal modes of small oscillations around 0o- Writing <p = O + Vi the normal 
modes are solutions of 

-^ + (V(x)+m 2 )r ] = uJ 2 r ] (2.3) 

with V(x) = U"((j)o(x)) — m 2 = §m 2 (0Q(:r) — 1). The quantum change in energy in 
going from the trivial vacuum to 4>o{x) is then 

£[M = \(H u - E "o) + £ct = AE + £ ct (2.4) 

where u are the free solutions (with (pl(x) = 1), and £ ct is the contribution from the 
counterterms. 

Note that it is possible that some of the values of oo 2 will be negative, so that 
in these directions our stationary point is a local maximum rather than a minimum 
of the action. These solutions will add an imaginary part to the energy, which we 
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can interpret via analytic continuation as giving the decay rates through the unstable 
modes [Q. If there is a direction in field space in which small oscillations lower the 
energy, we should be able to keep going in that direction and arrive at a lower min- 
imum. Later, we will see explicitly that the appearance of unstable modes coincides 
with the existence of a second solution to the </>o equation with the same J and lower 
energy. 

We would like to rewrite the sum over zero-point energies as an integral over 
phase space of the product of the energy and the density of states. We can then 
break this integral into a sum over bound states and an integral over a continuum, 
representing the latter in terms of phase shifts. In order to do so, however, we must 
review the peculiarities of Levinson's theorem in one dimension. For more details on 
these results see and 0. For a symmetric V(x), we can divide the continuum states 
into symmetric and antisymmetric channels, and then calculate the phase shift as a 
function of k separately for each channel, where uo 2 = k 2 + m 2 . The antisymmetric 
channel is completely equivalent to the I = case in three dimensions, so we have 

5a(0) = n A n (2.5) 

where Ua is the number of antisymmetric bound states. However, we must be careful 
in dealing with the special case of a state exactly at k = 0. In this case the solution 
to eq. ( |2.3|) with k = goes asymptotically to a constant as x — > oo, as opposed to 
the generic case where the k = solution goes to a constant plus linear terms in x. 
Just as in the I = case in three dimensions, this state contributes | to n&. We will 
refer to such states as "half-bound states." 

In the symmetric channel, Levinson's theorem becomes 

S S (0) = n s n-^ (2.6) 

where a bound state at k = contributes \ to ns, just as in the antisymmetric case. 
We can see the importance of getting the half-bound states right by looking at the 
free case: the phase shift is zero everywhere, and the right-hand side is zero because 
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the free case has a half-bound state (the wavef unction ip = constant). The situation 
is equally subtle for reflectionless potentials, all of which have half-bound states and 
Ss(0) + <5a(0) equal to an integer times n. 

We are now ready to rewrite the change in the zero-point energies in terms of 
phase shifts. Letting Ej be the bound state energies (again with k = bound states 
contributing with a |), p{k) be the density of states and po(k) be the free density of 
states, we have 



5?^-T + f^"s('*<" + «"> (27) 



where u(k) = y/k 2 + m 2 and we have used 



p(k) = p (k) + l^A(k) + 5 s (k)). (2.8) 



Note that the ^ term subtracts the contribution of the free half-bound state. 



Eq. (|2.7|) is divergent (the phase shifts fall as r for k — > oo), which is what we 



should expect since it includes the divergent contribution from the tadpole graph 
without the divergent contribution from the counterterms that cancels it. 

To avoid infrared problems later, we first use Levinson's theorem to compute the 
change in particle number, which is given by 

= El + jf f;£(fc(*> + *(*»-i (2-9) 

where again, half-bound states are counted with a | in the sum over j and the — \ 
comes from the contribution of the free half-bound state. Subtracting ^ times this 
equation from eq. ( |2.7|) , we have 

I /"OO fjh A 

AS = ^{E s - m) + j o -(u(k)- m )-(5 A (k) + 5 s (k)). (2.10) 
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Next, we and subtract the first Born approximation to eq. ( |2.10| ), which corre- 
sponds exactly to the tadpole graph. We must then add it back in using ordinary 
renormalized perturbation theory. However, in 1+1 dimensions we can adopt the 
simple renormalization condition that the counterterms cancel the tadpole graph and 
perform no additional finite renormalizations beyond this cancellation. With this 
choice, there is then nothing to add back in. 

The first Born approximation is given by 

(i) 1 f°° 

W = ~~T / V(x) cos 2 kxdx 



k Jo 

5^>(k) = --J V(x) sin 2 kxdx. (2.11) 

Notice that the sum of these two depends on V(x) only through the quantity (V) = 
Jq°° dxV(x), so we can indeed cancel the tadpole contribution with available countert- 
erms. 

Subtracting the first Born approximation, we have 



1 r°° dk d (V\ 

+ 1 ^(k)-m)-(5 A (k) + 5 s (k) + ^). (2.12) 



Since the Born approximation becomes exact at large k, the ^ term exactly cancels 
the leading | behavior of the phase shift at large k. As a result, this integral is 
completely finite and well-defined, since the integrand goes like -hs for k — > oo and 
goes to a constant at k = 0. We must find such a result, since in 1+1 dimensions 
eliminating the tadpole is sufficient to render the theory finite. 

In particular, we are free to integrate by parts, giving an expression that will be 
easier to deal with computationally, 

1 r°° dk k iV\ 

m = -J, sz^w+w + V'- (2 ' 13) 



Our use of Levinson's theorem to regularize AE in the infrared (replacing eq. (|2.7| ) 



by eq. ( |2.10p ) avoided a subtlety of the Born approximation in one dimension that 
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does not occur in higher dimensions: as k — > 0, the symmetric contribution from 
eq. ( j2.11j ) introduces a spurious infrared divergence, since it goes like |. (Each higher 
dimension adds a power of k near k = 0. We can see this in the Feynman diagram 
calculation, where these powers of k come from the measure.) 



2.2 Applications 



We can now use eq. (|2.13|) to calculate quantum energies for specific field configura- 



tions. Knowing that our model has a "kink" soliton solution <po(x) = tanh^rp, we 
will consider a family of field configurations that continuously interpolate between 
the trivial vacuum and a soliton/antisoliton pair, 



Tfl Tfl 

!>o(x, x ) = tanh — (x + x ) — tanh — (x — x ) — 1, (2-14) 



with 2xq measuring approximately the separation between the soliton and antisoliton. 
Unlike the kink, these configurations are not solutions of the equations of motion, 
except in the xq — > and xq — > oo limits. Thus we will need to to introduce a source 
that will vanish in these limits, and we must analyze the stability questions we raised 
earlier. 

For Xq very small, we simply have a small attractive perturbation from the trivial 
vacuum held in place by a small source, which we would not expect to introduce any 
instabilities. In terms of the scattering problem, for small xq, the potential is too 
weak to bind a state with a binding energy greater than m, which would give an 
imaginary eigenvalue. 

For Xq very large, we have a widely separated soliton/antisoliton pair. We know 
from translation invariance that a single soliton has a mode with u 2 = 0. Since this 
zero mode corresponds to a nodeless wavefunction, it is the lowest energy mode. The 
soliton/antisoliton pair has two translation modes that will mix, giving a symmetric 
eigenstate with a slightly lower energy and an antisymmetric eigenstate with a slightly 
higher energy. Thus we expect to find a single symmetric mode with uo 2 < 0, which 
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will contribute an imaginary part to £[0o]- The imaginary part gives the rate for 
our field configuration to decay through this mode toward the trivial vacuum. As 
xo — ► oo, we should find that the imaginary part goes to zero, and the real part goes 
to twice the energy of a single soliton, which we compute exactly below. 

According to this analysis, there should be a finite, nonzero value of Xq where 
the imaginary eigenvalue first appears. At this point, the field becomes unstable 
with respect to small perturbations in some direction in field space. Therefore, the 
energy must have a lower minimum that we can reach by moving in that direction. 
This configuration tp ( second stationary point of the action with the same 

source, which crosses (f>o at the value of xq where the imaginary part for the energy 
appears. (For smaller values of xq, this solution still exists but has higher energy.) 
In terms of the scattering problem, this crossing appears when the potential has a 
bound state with u 2 = 0. We can identify this state explicitly: since <po and ipo 
satisfy eq. (|2.2|) with the same J, as we approach the crossing, the wavefunction 
rj(x) = ipo(x, xo) — 4>o(x, xo) becomes a solution to eq. (|2.3|) with uj 2 = 0. 

We have carried out this computation numerically, and find results that agree 
with all of these expectations. To compute the antisymmetric phase shift, we will 
parameterize the wavefunction as 

v ( x ) = e~ ikx - e ikx e mx) (2.15) 

where (3(x) is an arbitrary complex function. Solving the differential equation for 
r)(x) subject to 7/(0) = 0, we find 

5 A (k) = -2Re/?(Jfe,0) ( 2 - 16 ) 
where (3(k, x) is at x = oo and satisfies 

- i/3" + 2k{3' + 2({3') 2 + ^V(x) = 0, (2.17) 
with prime denoting differentiation with respect to x. To obtain the symmetric phase 
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shift, we can follow the same derivation, but instead of imposing 77(0) = on the 
wavefunction, we instead impose r/(0) = 0, giving 



e -2i/3* k + 2/5. 
with (3 the same as above, so that 




Figure 2-1: The trajectory of £[4>o] in the complex energy plane, with £ in units of 
m, for £0 increasing from to — in steps of ^jp. Arrows indicate the flow along the 
real axis. 



Fig. |2- 1| shows the trajectory of £[(f>o] in the complex plane as a parametric function 
of xq, starting from the origin at Xq = 0. When xq ~ a single imaginary 

eigenvalue appears and £[<f>o] leaves the real axis. For xq large, £[4>o] approaches 
2(j^ — twice the standard result for a single soliton. The actual trajectory 
of £[4>o] has little significance, since it depends in detail on the functional form of 
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<p(x, Xq). However, the general features — beginning at the origin, moving up the real 
axis, out into the complex plane, and finally asymptotically to the real two-solution 
value — are characteristics of any <po that begins at the vacuum and ends at a well 
separated kink and antikink. 

Fig. |2-2| shows 4>q(x,Xq) and the second solution to the equations of motion with 
the same J, ip (x,x ). ipo goes to the trivial vacuum as x — > oo, and becomes a 
widely separated soliton/antisoliton pair as x — > 0, crossing O at Xo ~ 0.75. For x 
below this value, 0o has lower classical energy, while above this point, i/jq has lower 
classical energy, and this crossing appears precisely where the imaginary eigenvalue 
appears in the small oscillations spectrum. 




Figure 2-2: Solutions to the equations of motion with a source given by eq. (|2.2|) , 
as functions of x and xq in units of The shaded graph is 00(^,^0), which is 
guaranteed to be a solution by the construction of J. The unshaded graph gives the 
second solution ip Q (x,x ). (A finer mesh is used for this graph between x = and 
x ° = m or d er to illustrate its behavior in this region.) 

Through this continuous deformation from the trivial vacuum, we have arrived at a 
widely separated soliton/antisoliton pair, which we can now separate into independent 
configurations with nontrivial topology. These configurations are exactly soluble, so 
we will be able to study them analytically. However, our method does not rely on 
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having an analytic solution, so we could numerically calculate the energy of a generic 
field configuration with nontrivial topology using the same techniques. 



2.3 Analytic results 

In one dimensional quantum mechanics, potentials of the form 

V t (x) = -^m 2 sech 2 (^) (2.20) 

with £ an integer are exactly soluble and reflectionless. Their properties are sum- 
marized in Appendix A. The single soliton solution in our model, 4>o(x) = tanh^, 
corresponds to V(x) = — |m 2 sech 2 ^, the £ = 2 case of this family. (The sine-Gordon 
soliton corresponds to the £ = 1 case.) For a reflectionless potential, 5s(k) = 5\(k); to 
reconcile this equality with Levinson's theorem in the symmetric and antisymmetric 
channels, Vi(x) must have a half-bound state. (We saw this behavior already for the 
£ = case, the free particle.) We also note that although 5s = 5a, 5 s ^ 5 A , so the 
renormalized contributions from the symmetric and antisymmetric channels are not 
the same. In addition, the bound state contributions will also be unequal. 
For our case, the exact result for the phase shift is 

5s(A:) = 5AW = -tan- 1 -^^, (2.21) 

where the branch of the arctangent is chosen so that the phase shift is continuous 
and goes to zero for k — > oo. The Born approximation is 

6<»(k) + $\k) = *£. (2.22) 

There are three bound states: a translation mode with E — 0, a state with 
E = in ^-, and a half-bound state with E — m. Using eq. ( [2.12| ) or eq. (|2.13| ) we find 
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in agreement with [11 , [10], and the mode number cutoff method in M. 

We can calculate the energy of the soliton of the sine-Gordon model using the 

— I) 2 potential is replaced by 



same methods. In this model, the ^(0 2 



U(4>) = m 2 (cos 0—1) 



(2.24) 



which has soliton (anti-soliton) solutions 



4tan 1 (exp(=)=m(x — xq))). 



(2.25) 



The phase shift is 



8s(k) = 5 a (k) = tan 



-i 



rn 



(2.26) 



with Born approximation 



4 1) w+4 1, w = ^. 



(2.27) 



The bound states are just the translation mode at E = and the half-bound state 
at E — m, giving 

£[<h] = (2.28) 

7T 



again agreeing with the established results. 
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Chapter 3 



Fermions in one dimension 



We next extend our methods by adding fermions. We will see that our techniques 
extend cleanly and unambiguously to this case, and deal elegantly with the subtleties 
of boundary conditions. 

3.1 Formalism 

We consider a Majorana fermion \1/ interacting with a scalar background field 0, with 
the classical Lagrangian density 



where is the Lagrangian density for the background field, which we will take 
to be the same as in the last section, with the same soliton solutions. We note 
that choosing the Lagrangian density in this way causes the bosonic and fermionic 
degrees of freedom to be related by supersymmetry for a background field that is a 
solution to the equations of motion (such as a single soliton or an infinitely separated 
soliton/antisoliton pair). In this section, we will use this property only because we 
will find it instructive to compare the bosonic and fermionic small oscillations spectra. 
In the next section, we will consider the full supersymmetric model. 

The one-loop corrections to the energy due to fermionic fluctuations will be 




(3.1) 
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given by the appropriately renormalized sum of the zero-point energies, of 
the fermionic small oscillations. The spectrum of fermionic small oscillations in a 
background (f>o(x) is given by the Dirac equation, 

7° + V F (x)^j M%) = u£Mx) (3-2) 

where Vf(x) = m0 o (^) is the fermionic potential and k = ±V 'u 2 — m 2 is the momen- 
tum labeling the scattering states. 

We will choose the convention 7 = <r 2 and 7 1 = icr 3 so that the Majorana con- 
dition becomes simply = We note that since our Lagrangian is CP invariant, 
all of our results for the spectrum of a Majorana fermion can be extended to a Dirac 
fermion simply by doubling. 

Again, we will need to find the phase shifts. We can solve for the phase shifts 
of any fermionic potential Vp{x) that satisfies Vp(x) = Vp(—x) and Vp{x) -*mas 
x — > ±00 by generalizing the methods of the previous section. This form will be 
useful for considering our example of a sequence of background field configurations 
that continuously interpolates between the trivial vacuum and a widely separated 
soliton-antisoliton pair. In the limit of infinite separation, the phase shift for the 
pair goes to twice the result for a single soliton. For comparison, we also do the 
computation for a single soliton directly below. 

We define the parity operator acting on fermionic states as P = II7 , where II 
sends x — > —x. P commutes with the Hamiltonian, so parity is a good quantum 
number. As a result, we can separate the small oscillations into positive and negative 
parity channels, now restricted to k > in the continuum. 

We parameterize the fermion solutions by 

( e if(x) \ / e -iv(x)* \ 

e ikx and Wa;) = e~ ikx (3.3) 

where 9 = tan -1 — and v and £ are complex functions of x. We then find the phase 
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shift in each channel S ± (k) by solving eq. ( |3.2|) subject to the boundary conditions 

^(O)ocQJ (3.4) 

giving the phase shift as 

^±(0) = eT X2 (0) ± e^^e-f Xl (0). (3.5) 

Our boundary conditions assure that the wavefunctions ^ are eigenstates of the 
parity operator with eigenvalues ±1. We obtain 

5+(k) = -Re i/(0) + °- + log ^ ~ ' 



2 2i 1 — Y* 
5-(k) = -Rev^ + ^ + itogi^ (3.6) 

where F = ^(Vf(0) — ik + zi/(0)*) and z/(x) satisfies 

- zi/(x)" + v'(xf + 2Jfei/(x) + VF(a;) 2 - V F (x)' - m 2 = (3.7) 

with the boundary condition that v{x) and z/'(x) vanish at infinity. The total phase 
shift is given by summing the phase shifts in each channel, 5f(/c) = 5 + (k) + S~(k). 

Again, once we know the phase shifts, Levinson's theorem tells us how many 
bound states there will be. It works exactly the same way as in the bosonic case [0: 
In the odd parity channel, the number of bound states n_ is given by 

,T(0) = 7m_ (3.8) 

while in the even parity channel the number of bound states n + is given by 

5 + (0) = 7r(n + - \). (3.9) 

One can derive this result either by the same Jost function methods used in the boson 
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case, or by observing that at small k, the nonrelativistic approximation becomes valid 
so the bosonic results carry over directly. As in the boson case, for a particular 
potential there may exist a k = state in either of the two channels whose Dirac 
wavefunction goes to a constant spinor as x — > ±00. (Generically, for k = the 
components of the Dirac wavefunction go to straight lines as x — > ±00, but not lines 
with zero slope.) Just as in the bosonic case, such threshold states (which we will 
again call "half-bound states") should be counted with a factor of \ in Levinson's 
theorem. 

Given the phase shifts and bound state energies, we can calculate the one-loop 
fermionic correction to the energy. We continue to work in the simple renormalization 
scheme in which we add a counterterm proportional to <p 2 — 1, and perform no further 
renormalizations. The counterterm is fixed by requiring that the tadpole graph cancel. 
As in the bosonic case, we use the density of states 

pW=PoW + i^® (3.10) 

to write 

, u , l r f°° dk d5 F {k) 

£[M = -^i-J Q 2^-a- + £ >* ( 3 - n ) 

where u = \Jk 2 + m? and the sum over j counts bound states with appropriate factors 
of I as discussed above. Next, we use eq. (|3~lf) and eq. ( |3.9|) to rewrite this expression 

as 

„ r , ■, lr-i, \ f°° dk , ^dd F (k) 

In terms of the shifted field — 1, the (divergent) contribution from the tadpole 
graph is given by the leading Born approximation to 8p. A corresponding divergence 
related by <fi — > —<p symmetry appears in the second-order diagram. To subtract 
the tadpole graph, we subtract the leading Born approximation in Vp — m. In order 
to maintain the — > —0 symmetry, we also subtract the corresponding piece of 
the second-order diagram by subtracting the part of the second Born approximation 
related by the symmetry. Since we have chosen the counterterm to exactly cancel 
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these subtractions, there is nothing to add back in. Thus we have 

= -\ 5>i - m ) - r tr {u ~ m) Tk {5F{k) ~ 6?{k)) (3 - 13) 

where 

5 {1) (k) = ~f (Vf(x) 2 - V F (x)' - m 2 ) dx (3.14) 

which can also be obtained numerically by iterating eq. ( |3.7|) . For the soliton/antisoliton 
pair, the contribution from V'(x) vanishes since it is a total derivative, so that again 
our subtraction is indeed proportional to 2 — 1. 



3.2 Applications 

We continue to study a sequence of background fields labeled by a parameter xq that 
continuously interpolates from the trivial vacuum 4>(x) = 1 at xq = to a widely 
separated soliton/antisoliton pair at x — > oo, 

\ i mix — Xq) , mix + Xq) .„ , „. 

<j)(x, x ) = 1 + tanh —— — ^-tanh^— — ^. (3.15) 



Fig. |3-1| shows the fermionic bound state energies as a function of Xq. In the limit 
Xo — >• oo two bound states approach energy m^-. These are simply the odd and even 
parity combinations of the single soliton bound state at m^. The third (positive 
parity) bound state approaches uo = where the single soliton also has a bound state, 
but there is only one such state (with u > 0). Thus for a single soliton we must thus 
count the zero mode with a factor of |. We will see this result analytically below. 

For any finite xq, Levinson's theorem holds without subtleties; there are no states 
that require factors of |. Thus for a large but finite x , we find 

f<o) = I 

<T(0) = 7T (3.16) 
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Figure 3-1: Fermion bound state energies as a function of mx in units of m. Even 
parity states are denoted with <$> and odd parity states with +. 

consistent with having two positive parity and one negative parity bound states (see 



Fig. |3-1|) . In the limit Xq — > oo, an even parity "half-bound" threshold state enters 
the spectrum at uj = m just as in the bosonic case. Also, in this limit, the lowest 
(positive parity) mode approaches u = 0, and is only counted described 
above. Finally, a negative parity mode enters the spectrum from below, also to be 
weighted by ~. 

Thus the net change is to add half of a negative parity state, which via Levinson's 
theorem requires the phase shift <5~(0) to jump from n to ^ as Xq — ► oo. This jump 
occurs by the same continuous but nonuniform process as in all cases where a new 
state gets bound, which is illustrated in Fig. |3-2|. Just as in the bosonic case, in the 



limit of infinite separation the potential we have chosen becomes reflectionless, which 
requires 5 + {k) = 5~(k). 

To contrast the behavior of the zero modes, Fig. |3-3| shows u 2 for the bound states 
of the bosonic small oscillations. Because we have chosen the bosonic potential to be 
consistent with supersymmetry, the bosonic and fermionic spectra are related. Again 
as Xq — > oo the bound state energies approach those of the single soliton, and the 
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2.8 1 1 1 1 ' — 

0.05 0.1 0.15 0.2 

k/m 



Figure 3-2: Negative-parity phase shifts as functions of k/m for x = 2.0, 3.0, 4.0, 
5.0, 6.0, and 8.0. For any finite separation, the phase shift is equal to 7r at k = 0, but 
as xq gets larger, the phase shift ascends more and more steeply to 




_0.4 I i i i i i i i i i 

0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

mxo 

Figure 3-3: Boson bound state squared energies as a function of mxo in units of m. 
Symmetric states are denoted with <J> and antisymmetric states with +. 
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wavefunctions are formed from the odd and even combinations of the wavefunctions 
for the single soliton. However, we see that in the boson case both the mode at m^- 
and the mode at uj = are doubled, so there is no factor of | in counting the bosonic 
zero modes. 

0.4 



AH 




Figure 3-4: One-loop fermionic correction to the energy as a function of Xq in units 
of m. 



Fig. 3j4 shows the values of £[0o] we obtain from eq. ( 3.13|) as a function of Xq. 
In the Xo — > oo limit, we can also do the calculation analytically. We will expand on 
this technique in our discussion of the supersymmetric case below. We find a phase 
shift 

(3.17) 



S F (k) = 4 tan" 1 — + 2 tan" 1 — 
2k k 



with Born approximation 



(3.18) 



and contributions from two bound states at u — m^- and one at to = m. We thus 



find 



£\ 



2m 



1 



which agrees with the numerical calculation. 



(3.19) 
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Chapter 4 

The supersymmetric model in one 
dimension 



Having analyzed bosons and fermions separately, we now combine these results and 
consider the full supersymmetric mode in one dimension. Here we will work analyt- 
ically using a general potential, with the 4 and sine-Gordon theories providing two 
concrete examples. The supersymmetric model also introduces another quantity, the 
central charge, which we must also compute in order to check the consistency of our 
results for the energy. A review of the literature shows a wide variety of conflicting 
results]^, [T2|, |T^, |T3[|. Our techniques for dealing with the regularization, renor- 



malization and calculation of one-loop corrections to the energies of classical field 
configurations will enable us to take a fresh look at this problem. We will obtain an 
unambiguous result for the energy and the central charge, and in the process we will 
see how our techniques are related to the expansion of the quantum field in creation 
and annihilation operators. 

4.1 Formalism 

We begin with the classical supersymmetric Lagrangian density in 1+1 dimensions 

2 

C = J\ ((^X^) - W) 2 + - u'{4>)^) (4.1) 
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where is a real scalar, \1> is a Majorana fermion, and U((f>) = y(0 2 — 1) for the 4 
model and U(4>) = 2msin(0/2) for the sine-Gordon model. These models support 
classical soliton and antisoliton solutions, which are the solutions to 

<f>' (x) = T U(<j> ) (4.2) 

for the soliton and antisoliton respectively. In the 4 model, the soliton solution is 
the "kink" that we have already discussed, 

771 T 

0^ k (x)=tanh— , (4.3) 
while the soliton in the sine-Gordon model is given by 

</) S G {x) =4tan- 1 e- mx , (4.4) 

and in both cases the antisoliton is obtained from the soliton by sending x — > —x. 
The eigenvalue equations for the bosonic and fermionic small oscillations are 

' <Jl +U'(<p ) 2 + U(fo)U"^ ))r ]k (x) = (cu*) 2 Vk (x) (4.5) 



dx 2 



7 °(-V^ + ^U))^(z) = w^ fc (x). (4.6) 



The bosonic potentials are given by 

£/'(0o) 2 + Wo)t/"(0o) -m 2 = - ^±^j m 2 sech 2 ^ = V e (x) (4.7) 

with £ = 1 for the sine-Gordon soliton and £ = 2 for the kink. Repeating our results 
from above, we have phase shifts 

5^(k) = 5 e=2 (k) = 2 tan" 1 '' 3mk ^ 



2k 2 — m 2 
m 



S s B G (k) = S e=1 (k) = 2tan~ 1 ^ (4.8) 
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and bound states at u = for the kink ui = ;y y 2: , and at u = for the sine-Gordon 
soliton. As discussed above, both potentials have half-bound states at threshold, a 
necessary consequence of their reflect ionless property. 

As above, we avoid subtleties of boundary conditions by computing the fermionic 
phase shifts for a widely separated soliton and antisoliton solution, with the antisoliton 
on the left so that U'(4>) — > m as x — > ±00. We will use the second-order equation 
obtained by squaring eq. (|4.6|), 



tl) k {x) = k 2 Mx) (4.9) 



for the soliton and 



"Sr + Vt(x) 

~& + Vt(x) 



-& + Vt(x) 

~£ + V e (x) 



■4> k {x) = k 2 M%) (4-10) 



for the antisoliton, where Ve(x) = ^V^-i (f ). An incident wave far to the left is given 
by 

* W " (4 ' U) 

with 6 = tan -1 — . To start with, we will restrict to the case of a reflectionless 
potential, such as the kink or sine-Gordon soliton. It thus scatters without reflection 
through the antisoliton becoming 

where 6i(k) is the phase shift for the bosonic potential V^{x) and S^(k) is the phase 
shift for the bosonic potential Ve(x). It then scatters without reflection through the 
soliton giving 

M*) = * ikx \ . t ,t ,x ,« I- ( 4 - 13 ) 
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By rescaling the results from Appendix A, we easily obtain 



~ 771, 

S t (k) = S t (k) - 2 taiT 1 — (4.14) 



so that for the soliton/antisoliton pair, 



171 

S F (k) = 5 B (k) - 2 tan" 1 — (4.15) 

k 



and thus for a single soliton 



Til 

S F (k) =5 B (k) -tan" 1 — (4.16) 

K 



in both models. This result has been also been obtained in |T3| and M. Through 



this analysis, we see that the deficit between the boson and fermion phase shifts is 
necessary so that eq. ( 4. 12|) correctly solves the Dirac equation in the region where 



Vf(x) = —m. This result also agrees with the methods of the previous section. 

These results generalize to any supersymmetric potential U(4>) that supports a 

soliton 0o with 4>o(x) = —(f>o(—x). We can still consider eq. (|4.9| ), with Vt(x) and 
Vt(x) replaced by 

V{x) = U'ifa) 2 + UMU'Xfa) - m 2 (4.17) 

and 

V(x) = U'(<p ) 2 - C/(0q)[/"(0 o ) - m 2 . (4.18) 

These are symmetric, though now not necessarily reflectionless, bosonic potentials. 
We decompose their solutions into symmetric and antisymmetric channels. For x — > 
±cx3, these solutions are given in terms of phase shifts as 

r)l(x) = cos(kx ± 5 s (k)) Vk( x ) = sin(A;x ± 5 A (k)) 

rj%(x) = i cos(kx ± 5 s (k)) fj£(x) = -i sin(A;x ± 6 A (k)) (4.19) 

where the arbitrary factors of ±i are inserted for convenience later. For all x these 
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wavefunctions are related by 

^vl(x) = i {J- + C/'(0 O )) rt{x) u kV A (x) = i(J^ - f/'(0 o )) f k (x) 

u k rii(x) = i ( ^ - U\</>o)) fj£(x) 



w k fj£(x)=i^ + U'((f> ) y j 

(4.20) 

so that the solutions to the Dirac equation are 



dx 



ri(x)=( Vk A ) and ^(x)=h) (4.21) 

with positive and negative parity respectively. The phase relation between the upper 
and lower components of these wavefunctions must be different as x — > ±00 since the 
mass term has opposite signs in these two limits. 
Putting this all together gives, as x — > ±00, 

cos(kx±5 s (k)) = — ( 4- ~ V F (x) ) sin(kx±5 A (k)) 

LO k \dX 



— — (k cos(kx ± 5 A (k)) msm(kx ± 5 A (k 
u k v 



= t sin(kx ± 5 A (k) t0) = cos(kx ± 5 A (k) T ±^-) (4.22) 



and thus 



and similarly 



5 s (k) = 5 A (k)+ tan' 1 ^ (4.23) 

K 



S A (k) = 5 s (k)+tan~ 1 ^. (4.24) 

K 



The fermion phase shift in each channel is given by the average of the bosonic phase 
shifts of the two components 

5 + (k) = -(5 s (k)+5 A (k)) = 5 s (k)--tan- 1 ^- 

2 2 rC 

S~(k) = \{5 A {k) + ~5 s {k)) = 5 A {k)- l -tan^ 

(4.25) 
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so that 



5 F (k) = 5 + (k) + 6-(k) = 5 s (k) + 5 A (k) - tan" 1 ? (4.26) 

k 



and we obtain the same result, eq. (|4.16 ), as we found in the reflectionless case 



Since 5p(0) ^ Sb(0), Levinson's theorem requires that the spectrum of fermionic 
and boson bound states differ. The difference is that, although there is a fermionic 
bound state for every bosonic bound state, the mode at uj = only counts as | for the 
fermions. (The fermionic states at threshold also count as ~, the same as in the boson 
case.) We have seen this difference in the previous section by the tracking the bound 
state energies as we interpolated between the trivial vacuum and a soliton/antisoliton 
pair. We can also check this result analytically by observing that the residue of the 
pole at k — im in the reflection coefficient Tp is half the residue of the pole at k — im 
in Tb because of eq. (|4.16|) . As a final check, we imagine doubling the spectrum 
by turning into a complex scalar and \I/ into a Dirac fermion. Then in a soliton 
background <fi would have two zero-energy bound states, one involving its real part and 
one involving its imaginary part. However, \1/ would have only a single zero-energy 
bound state, with wavefunction given by 




1>(x) = (4.27) 



with Vp(x) = U'(<fto). The corresponding solution with only an lower component is 
not normalizable; for an antisoliton, we would find the same situation with upper and 
lower components reversed. Thus when we reduce to a Majorana fermion, we count 
this state as a half. 

We note that the fermionic phase shift and bound state spectrum are simply 
given by the average of the results we would obtain for the two bosonic potentials 
Ve(x) and Vg(x). We also note that just as the bosonic zero mode arises because the 
soliton breaks translation invariance, the fermionic zero mode arises as a consequence 
of broken supersymmetry invariance (which we can think of as breaking translation 
invariance in a fermionic direction in superspace). For a soliton solution, only the 
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super symmetry generator Q_ is broken, while Q + is left unbroken (the situation is 
reversed for an antisoliton) . Thus since the supersymmetry is only half broken, it is 
not surprising that the corresponding zero mode counts only as a half. In both cases, 
acting with the broken generator on the soliton solution gives the corresponding zero 
mode. 



4.2 Applications 

To compute the one-loop correction to the energy, we now follow the same method 
as in the previous sections and sum the quantity |u; over bosonic and fermionic 
states, with the fermions entering with a minus sign as usual. We will discuss the 
case of an isolated soliton, and see that results agree with the widely separated soli- 
ton/antisoliton pair considered in the last section. 

Thus our formal expression for the energy correction is 



where the states at threshold and the fermion bound state at u = are weighted by 
| as discussed above. The free density of states has cancelled between bosons and 
fermions, as required by supersymmetry. Again, to avoid infrared problems later, we 
use Levinson's theorem 

5(0) = n(n B - ~) (4.29) 

to rewrite eq. ( |4.28| ) as 

* r-r 1 x—w r x lv^/ p x f°° dk , . / (15b d5p \ , , „^ 

= 5 XX - "0 " 2 XX - H + 1 S <« - "0 (-a— J <«°) 



where the \ in eq. ( |4.29| ) has cancelled between bosons and fermions. 



The continuum integral in eq. ( |4.30| ) is still logarithmically divergent at large 
k, as we should expect since we have not yet included the contribution from the 
counterterm. As above, we can isolate this divergence in the contributions from the 
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low order Born approximations to the phase shifts 5b and 8p. We then identify these 
contributions with specific Feynman graphs, subtract the Born approximations, and 
add back in the associated graphs. For the boson, the divergence comes from the first 
Born approximation, which corresponds exactly to the tadpole graphs with a bosonic 
loop. For the fermion, the source of the divergence is more complicated: we subtract 
the first Born approximation to the fermionic phase shift and the piece of the second 
Born approximation that is related to it by the spontaneous symmetry breaking of 
<f). This subtraction corresponds exactly to subtracting the tadpole graph with a 
fermionic loop and the part of the graph with two external bosons and a fermionic 
loop that is related to the tadpole graph by spontaneous symmetry breaking (the 
rest of the two-point function is then finite). For both the boson and fermion, this 
subtraction amounts to simply subtracting the term proportional to | that cancels 
the leading \ behavior of the phase shift at large k. We can identify the coefficient of 
these | terms with the coefficients of the logarithmic divergences in the corresponding 
diagrams. As a result, by computing the divergences in the bosonic and fermionic 
diagrams, we obtain a check on eq. ( |4.16|) , to leading order in r for k large. 



Of course we must add back all that we have subtracted, together with the con- 
tribution from the counterterm. To do so we must consider renormalization. We will 
continue to use our simple renormalization scheme, in which we introduce only the 
subtraction 

C^C- CU"(<j))U{<f)) - CU"'{(j))^ (4.31) 

which is equivalent to 

U{(j>) -> U(<j>) + \cu"{<j>) (4.32) 

and thus preserves supersymmetry. We fix the coefficient C by requiring that the 
boson tadpole (which includes contributions from both boson and fermion loops as 
we have described above) vanish. In this scheme, the counterterm completely cancels 
the terms we have subtracted from eq. ( |4.30| ), so there is nothing to add back in. In 
the sine-Gordon theory, this scheme also makes the physical mass of the boson equal 
to m, while in the 4 theory, there is a one-loop correction to the physical mass of 
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the boson from the diagram with two three-boson vertices, giving a physical mass 
of m — ||. For us it is more important to guarantee that the tadpole graphs 

vanish, assuring us that we have chosen the correct vacuum for the theory, than to 
have the physical mass equal to the Lagrangian parameter m; for the sine-Gordon 
case we happen to be able to do both at once. Furthermore, such renormalization 
conditions can be applied uniformly to arbitrary U{4>). 

Thus the effect of regularization and renormalization in our renormalization scheme 
is to subtract 

5^(k) = 5^(k)-5?(k) = j (4.33) 
from the difference of the boson and fermion phase shifts, giving 

a# = \ HK B ~ m )-\ HK F - m ) 

j 3 

r°°dk, Jd8 B dd F d5M\ 

+ Jo ^- m) \-dk-lk-^k) 

m r°° dk , , d ( i m m\ m 



"4 + 



f°° dk . .d ( ,171 m\ m , Jn ^ 

/ — (u - m)— - tan" 1 — - — = (4.34) 

Jo 2n J dk\ k kj 2vr v ; 



for both the kink and sine-Gordon soliton. This result agrees with H and [14], and 



disagrees with |]I2"|| , p[, |T3]| , and As pointed out in 0, in the case of the sine- 
Gordon soliton, it also agrees with the result obtained from the Yang-Baxter equation 
assuming the factorization of the S- matrix ||17|| . 

We note that in the end this result depended only on eq. ( |4.16| ) and its impli- 



cations for Levinson's theorem. Thus, since eq. ( [4.161) holds in general for antisym- 



metric soliton solutions, eq. ( |4.34j ) gives the one-loop correction to the energy in our 
renormalization scheme of any supersymmetric soliton that is antisymmetric under 
reflection. 



4.3 Supersymmetry algebra and the central charge 

Our second application of the apparatus we have developed is to compute the one-loop 
quantum correction to the central charge in the presence of the kink or sine-Gordon 
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soliton. 

First we summarize the supersymmetry algebra. We define 



Q ± = { T 11 J ^- / {p(j> + iU)^ dx = ^- f (n*± + (</>' ± rfx (4.35) 

2 \ J X J 

where = 1=F j? ^ and Q± = 1=F ^ 71 Q. Using the canonical equal-time (anti) commutation 
relations, we have 

TV? 

— {i*±{x), *±(y)} = i5(a: - y) 

2 

— [0(x), n(y)]=i5(x-y) (4.36) 

where II = is the momentum conjugate to <fi and all other (anti) commutators vanish. 
The supersymmetry algebra is 



{Q±, Q±} = 2H±2Z {Q + ,Q_} = 2P, (4.37) 
where H, P, and Z are given classically by 

2 

H = ^- J (IT 2 + (0') 2 + U 2 + + + 2zf/'^_^ + ) 

P = ^|(n0' + ^(*^V + *-*'-)) 

9 ~" 

777 /" 

Z = —JfiUdx. (4.38) 

It is easy to check that H is the same Hamiltonian as would be determined canonically 
from eq. ( |4.1|) . 

At the classical level, using eq. (|4.2| ). 



2 2 

Hcl = ! kl ( '° (X)2 + ^ (0o)2 ) dX = T Tl U (MxMW = TZ ch (4.39) 

for the soliton and antisoliton respectively. 

The hermiticity of Q± gives the BPS bound on the expectation values of H and 
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Z in any quantum state: 

(H)> \(Z)\. (4.40) 

Classically, the values of H and \Z\ are equal so this bound is saturated. We have 
found a negative correction to H at one-loop, so if there is no correction to Z, eq. ( [4.40[ ) 
will be violated. 

To unambiguously compute the corrections to the central charge for a soliton, 
it is easier to consider corrections to Q\ = H + Z, which is zero classically (for the 
antisoliton we should consider Q 2 _ = H+Z). One reason to consider Q 2 + rather than H 
and Z separately is that this quantity is finite and independent of the renormalization 
scheme. Using eq. (|4.31p and eq. (f4.32p we see explicitly that the contribution from 



the counterterm cancels: 

AH ct = C J U"(<t> )U((/> ) dx = -C J f/"(0 o )0[, dx = -AZ ct (4.41) 

(and we only need consider the tree-level contribution since the counterterm coefficient 
C is already order A ). 

Next we expand <f>(x) = 4>o(x) +r](x), where the soliton solution (fio is an ordinary 
real function of x. Neglecting terms of order rf and higher (which give higher-loop 
corrections), we obtain 



2 

m 



(h + zu = ^/(n 2 + 



(J^ + tffafjV+J dx } (4.42) 

where ()^ denotes expectation value in the classical soliton background. 

To evaluate this expression, we decompose the fields rj and $ using creation and 
annihilation operators for the small oscillations around <pQ. The small oscillation 
modes will be given in terms of the eigenmodes of the bosonic potentials Vg(x) and 
Ve(x) = -kVe-i(j)- For any mode r/fc(x) of Vg{x) with nonzero energy Uk = y/k 2 + m 2 , 
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there is a mode fjk(x) of Vf(x) with the same energy, related by 



UkVk(x) = « 7 + U'((j) ) j r/fe 



WA77 fc (z) = if ^ - U'{<f) ) ) fj h . (4.43) 



We use these wavefunctions to obtain 



" w = V^/7i^( <, * (l)e_ '"* , + aI "' (x)e " 

+\/~V^=o( :r ) a ^=o 

V_?Ti 

+\/^ w =o(0)6 w=0 (4.44) 

V 771 

where 77_fc(:r) = i]l(x), the creation and annihilation operators obey 

with all other (anti) commutators vanishing, and 

(r/fc (x) \ 
(4.46) 
%0)/ 

We note that the integral over k also includes discrete contributions from the bound 
states (which correspond to imaginary values of k). These are understood to give 
discrete contributions to the results that follow (with Dirac delta functions replaced 
by Kronecker delta functions appropriately). However, we have explicitly indicated 



the contribution from the bound states at uo = following |T8 
We normalize the wavefunctions rjk such that 



dk 

—Vk(x)*r]k(y) = S(x - y) (4.47) 

Z7T 
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which implies 

7^Vk(x)*fj k (y) = 8(x-y). (4.48) 

With this normalization, the fields 77 and \I/ obey canonical commutation relations. 
Elementary algebra yields 




dx J —u k \fj k {x)\ 2 - J dx J —u k \r] k {x)\ 2 = (4.50) 

and the BPS bound remains saturated. (If we instead considered an antisoliton, 
we would find the same result for {Qt)g = (H — Z)z, with the roles of and 
ty- reversed.) Our result disagrees with |§ and [|15f| , which claim that there is no 
correction to the central charge at one loop in this renormalization scheme. We note 
that this result did not depend on any specific properties of U, so it holds for any 
super symmetric soliton satisfying eq. ( |4.2|) . 

The second line of eq. ( |4.50| ) is simply the unregulated fermionic contribution to 
the energy, and is explicitly equal to minus the average of the contributions from 
the bosonic potentials Vi and Vt, in agreement with what we found above. As a 
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final consistency check, we recalculate the full one-loop correction to the energy and 
central charge using our expansion in terms of quantum fields. For AH we obtain, 
again neglecting rf terms, 



AH = (H)4-H d 

= ^ / (n 2 + v [~^ 2 + f/'(0 o ) 2 + Wo)t/"(0o)) v 

= AH ct 

+ J dx J —u k \i] k (x)\ 2 - J dx J —u k (\rj k (x)\ 2 + \i] k (x)\ 2 ^j . (4.51) 

To relate this expression to our phase shift formalism, we consider the Green's function 
for the bosonic field 



G(x,y,t) = iT (rj(x,t)r](y,0)) 

^ +e-^r ]k (x)r ] t(y)e(-t)) (4.52) 



and its Fourier transform 



G(x, y, U ) = / G(x, y, t) e «* ={£■( \ t(x) f y) ) (4.53) 

J J Z7T \ Ld z — UJJ, — It I 



whose trace gives the density of states according to 

p B {uo) = \m — / G(x,x,uj) dx (4.54) 

7T J 

giving as a result 

p B (k) = - I dx\ Vk (x)\ 2 . (4.55) 

71 J 

Similarly for the fermions we find 

p F (k) = ^Jdx (\rj k (x)\ 2 + Mx)\ 2 ) ■ (4-56) 
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These results enable us to verify that eq. Q4.51 ) is in agreement with eq. ( 4.34 ). 



In the exact same way, we can calculate the correction to Z directly. We start from 



the classical expression for Z in eq. (|4.38|) and expand about the classical solution 
<P = 0o, giving 



AZ = (Z) 4> -Z t 



cl 

2 

m 



AZ ct + — 



j (u'ryrf - ifU"rj^ dx 

AZct + / ( { { 1 + u ' )v ) 2 - iv ' )2 ~ v2(u ' )2 ~ uu " v2 ) dx - {4 - 57) 



After substituting the expansions of eq. (|4.49|) we obtain 



AZ = AZ ct + J dx J ^-u k \r]k(x)\ 2 - J dx J ^-u k \r] k (x)\ 2 

3 3 

+ J ^-(uj k - m)^ ($i(k) - 8i(k) + 25<»(k)) 

= / — (uo-m — tan" 1 — - — = — = -AH. 4.58 

4 J 2vr v J dk\ k k J 2tt v 1 



Another work [20] has explained the difference between our result for the central 
charge and the result obtained in Ref. || and predecessors. The earlier works found 
AZ = based on an argument that involves direct manipulation of the Z operator. 
However, Ref. [^0] showed that manipulations of this type are only valid if the operator 
is augmented with an anomalous correction of order h. 

We have systematically avoided the questionable manipulations that would have 
led to AZ = by computing only matrix elements, where we have the phase shift 
formalism available to guide us. First we developed an unambiguous renormaliza- 
tion procedure for H based on physical quantities. Classical BPS saturation defines 
the operator Z at tree level. Once we have fixed a renormalization scheme in our 
computation of H, the expectation value of Z (or any other physical quantity) is 
determined. We then carried out the one-loop computation of Z exactly in parallel 
to the computation of H so that no new ambiguities could arise, and found that BPS 
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saturation is maintained at one loop, as confirmed in Ref. [20]. Our use of Levinson's 
theorem and the Born approximation also prevented the appearance of the spurious 
linear divergences found in Ref. [§]. 

Thus we consistently included the effects of the SVV anomaly in both H and 
Z, though of course we found only a particular matrix element instead of the full 
operator. However, by working in the continuum (without boundaries), we avoided 
the difficulties that Ref. |2(| faced in separating the unphysical contributions at the 
boundaries from the physical effects localized at the soliton. 
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Chapter 5 

Scalars in three dimensions 



Having fully analyzed models in one dimension, we now turn to three dimensions. 
We will continue to consider models with spherical symmetry, so that we can use a 
partial wave decomposition. Thus, in addition to an integral over k, we will also have 
a sum over £ to consider. On the other hand, the subtleties of the symmetric channel 
in one dimension will be absent. Theories in three dimensions will also have stronger 
divergences than we found in one dimension. Thus we will have to subtract more Born 
terms and add back in more corresponding diagrams. However, since all the theories 
we will consider are renormalizable, only a finite number of Born subtractions will 
be required. We will again begin with a scalar model, but our techniques will easily 
generalize to include fermions. 

5.1 Formalism 

We consider a renormalizable field theory with a real scalar field coupled to a 
charged scalar \1>. We take the classical potential V{4>) oc ((f) 2 — v 2 ) 2 , and ^ acquires 
a mass through spontaneous symmetry breaking. At the quantum level we put aside 
the self-couplings and consider only the effects of the <fi — \l/ interactions. We further 
restrict ourselves to 0(K) effects in the quantum theory, which correspond to one-loop 
diagrams. 
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Our model is defined by the classical action 



J d 4 x |i(9 M 0) 2 - ±{<f _ v 2y + g^dHx _ gy*^ 
+a{d^) 2 - b(<p 2 - v 2 ) - c(0 2 - v 2 ) 2 } , 



where we have separated out the three counterterms necessary for renormalization 
and written them in a convenient form. At one-loop order in these are the only 
counterterms required. 

We quantize around the classical vacuum <p = v and define h = <fi — v , so that 

S[h,V] = jd 4 x^(d ll h) 2 -^(h 2 + 2vh) 2 -g(h 2 + 2vh)^ 
+d^*d^ - M 2 i)*m 
+a{d^h) 2 - b(h 2 + 2hv) - cih 2 + 2hvf | 

where M = yjgv is the \& mass and m 2 = Xv 2 /3 is the h mass. 

The one-loop quantum effective action for h is obtained by integrating out \& 
to leading order in H. We are interested in time-independent field configurations 
h = h(x), for which the effective action yields an effective energy £[h] that has three 
parts: 

£[h] = £ cl [h] + £ ct [h] + £*[h] , (5.1) 
where £ c \[h] is the classical energy of h, 

£ d [h] = J d 3 x ji|V/f + ^- 2 {h 2 + 2vh) 2 ^ , (5.2) 

£ c t[h] is the counterterm contribution, 

S ct [h] = J d 3 x [a\Vh\ 2 + b{h 2 + 2hv) + c{h 2 + 2hv) 2 ) , (5.3) 

and £$[h] is the one-loop quantum contribution from \E'. £ ct [h] and £^[h] are divergent, 
but we will see explicitly that these divergences cancel for any configuration h(x). 
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2gvtt 



2vb 



(2gvf U(p 2 " 



g n 



b + Av c ap 
Figure 5-1: One- loop diagrams. 

We fix the counterterms by applying renormalization conditions in the perturba- 
tive sector of the theory. Having done so, we have defined the theory for all h(x). We 
choose the on-shell renormalization conditions 



Si = 0, S 2 (m 2 ) = 0, and 



dS 5 



dp 2 



0. 



(5.4) 



where Si and S 2 (p 2 ) are the one- and two-point functions arising only from the loop 
and counterterms as seen in Fig. |5-lj We denote the one-loop diagrams with one 
insertion by Q and with two insertions by U(p 2 ), and find 



Si = 2vgtt + 2vb , 
£ 2 (jo 2 ) = (2vg) 2 U(p 2 ) + gtt + b + (2v) 2 c + ap 2 . (5.5) 



Defining 



nV) = 



dU(p 2 
dp 2 



the renormalization conditions eq. (|5.4j ) then yield 



(5.6) 



a = -(2vgyW(m 2 ) , b = -gQ , c = g 2 (m 2 W(m 2 ) - U(m 2 )) 



(5.7) 



which we then substitute into the counterterm energy, eq. (5.3). 

Now we consider the calculation of £q, [h] . This energy is the sum over zero point 
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energies, \fioJ, of the modes of \1/ in the presence of h(x), 

= $>*[/i] (5.8) 

a 

where ouj~ are the positive square roots of the eigenvalues of the small oscillations 
Hamiltonian, 

' _v 2 + M 2 + g{h 2 + 2vh)) i) k = u 2 a ij k (5.9) 



The fact that \I> is complex accounts for the absence of | in eq. (|5.8|) . 

is highly divergent. However our model is renormalizable and therefore the 
counterterms fixed in the presence of the trivial h must cancel all diverg ences in £^j. 
Rather than attempt to regulate the divergent sum in eq. ( |5.8| ) directly, we study 
the density of states that defines the sum. We can isolate the terms that lead to 
divergences in and renormalize them using conventional methods. Thus our task 
is to generalize the construction of eq. ( |2.7|) to this case. 

For fixed h(x) the spectrum of H given in eq. (513) consists of a finite number 



(possibly zero) of normalizable bound states and a continuum beginning at M 2 , pa- 
rameterized by k, with E(k) = +y/k 2 + M 2 . Furthermore, H depends on h only 
through the combination 

X = h 2 + 2hv , (5.10) 

so we can consider to be a functional of x- We restrict ourselves to spherically 
symmetric h. Then 

£*M =E^' + E(2^+l) / dk Pe (k)E(k) (5.11) 

j £ 

where pe(k) is the density of states in k in the £ th partial wave and the Ej are the 
bound state energies. pe{k) is finite, but the sum over I and the integral over k are 
divergent. Furthermore 

pt(k)= P r c (k) + 1 -^± } (5.i2) 

where 5e(k) is the usual scattering phase shift for the £ th partial wave, and pf ee (k) is 
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the free (g = 0) density of states. This relationship between the density of states and 
the derivative of the phase shift is shown for example in ]TJJ. 

At the outset, we subtract p free (/c) from the density of states since we wish to 
compare £q,[x] to £#[0]. Viewing as the sum of one loop diagrams, we see that 

only the diagrams with one or two insertions of g\ are divergent. A diagram with n 
insertions corresponds to the n th term in the Born expansion, so all possible diver- 
gences can be eliminated by subtracting the first and second Born approximations 
from the phase shifts that determine the density of states. Standard methods allow 
us to construct the Born approximation for the phase shifts ||, which is a power 
series in the "potential" g\- 

We define the combination 

5 e (k) = S t (k) - 5f\k) - 8f\k) , (5.13) 

where 5 t (k) and Sf\k) are the first and second Born approximations to 5i(k). We 
then have 

*W = E^ + E(« + Djf «i^(*) + *n/^i«o) (5.i4) 

where 

X(p) = J d 3 x X (x)e- 1 ^ , (5.15) 

and likewise for h(p). Both h and x are real and depend only on q = \p\ for spherically 
symmetric h. We have subtracted out the order g and g 2 contributions by using 
Se(k) instead of 5e(k), and added them back in by using their explicit diagrammatic 
representation in terms of the divergent constant Q and the divergent function IT(p 2 ). 
We can now combine and £ ct and obtain a finite result: 

£* + £ ct = $>,+£(2£ + l) rdk-^^-E(k) + T 2 [h} (5.16) 
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where 



r a [fc] 



q 2 dq - 



9 M 



2 Jo 2tt 2 
' T ./) 2tt 2 



q 2 dq 



(n(-g 2 ) - n(m 2 ) + m 2 n'(m 2 )) x(g) 2 

4v 2 g 2 n , (-g 2 )/i(g) 21 



(5.17) 



II is log divergent, but both {II(— q 2 ) — Il(m 2 )} and II' are finite, so T 2 [h] is finite as 
well. 

Each term in the Born approximation to the phase shift goes to zero at k = 0, so 
by Levinson's theorem 5e(0) = Se{0) = Ttng where ri£ is the number of bound states 
with angular momentum £. As k — > oo, 5e(k) falls off like i 5^\k) falls off like |, 

(2) l 

and 5} (k) falls off like p. Since the Born approximation becomes exact at large fc, 
5^(/c) falls like p-. Thus we see that the first subtraction renders each integral over k 
convergent. The second subtraction makes the f-sum convergent. We are then free 
to integrate by parts in ( J5.16 ), obtaining 



1 rc 

e[h] = £a[h] + r 2 [h]-~ y E(u + i) / 

7T c JO 



dk Se(k) 



k 



E{k) 



+ J2(E J -M). {5.11 



In this expression we see that each bound state contributes its binding energy, Ej —M, 
so that the energy varies smoothly as we strengthen h and bind more states. 



5.2 Calculational methods 



In this Section we describe the method that allows us to construct £ [h] as a functional 
of h and search for stationary points. We now consider the calculation of each of the 
terms in eq. (|5.18 ) in turn. The classical contribution to the action is evaluated 
directly by substitution into eq. ( |5.2|) . I^f/i] of eq. (|5.17 ) is obtained from a Feynman 
diagram calculation, 



T 2 [h] 



(4tt) 



q 2 dq 
(2tt) 2 



log 



m x(l — x) 



M 2 — m 2 x(l — x) 



M 2 + q 2 x(l - x) 
M 2 — m 2 x(l — x) 

x(qfdx 
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q 2 dq 



(4vr) 2 Jo (2tt) 2 Jo 



x(l — x) 



M 2 — m 2 x(l — x) 



Av 2 q 2 h(qY 



dx. (5.19) 



The partial wave phase shifts and Born approximations are calculated as follows. 
The radial wave equation is 



< + 



£(£ + !) 



+ gx{r) 



Ui = k 2 u e , 



(5.20) 



where k 2 > 0, and %(r) — > as r — > oo. We introduce two linearly independent 
solutions to eq. (|5.20|) , (r) and uf\r), defined by 



uP(r) = e 2lMk ' r) rh ( l\kr) 
u {2 \r) = e- 2i ^ {k ' r) rhf\kr) 



(5.21) 



where hfp is the spherical Hankel function asymptotic to e bkr jr as r — > oo, hf\kr) = 
hP*{kr), and (3e(k,r) — > as r — > oo, so that u\ (r) — > e ifcr and (r) — > e~ lfcr as 
r — > oo. The scattering solution is then 



.(2) 



uAr) = u) (r) + e n ; w) (r 



2^(fc)„(i), 



(5.22) 



and obeys m^(0) = 0. Thus we obtain 



8t(k) = -2Rep e (k,0). 



(5.23) 



Furthermore, fit obeys a simple, non-linear differential equation obtained by substi- 
tuting vJjP into eq. ( |5.20|) , 



- - 2tkp e (kr)(3' e + 2((3' e ) 2 + -g X (r) = 



(5.24) 



where primes denote differentiation with respect to r, and 



xh^\x) 



(5.25) 
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is a simple rational function of x. 

We solve eq. ( |5.24j ) numerically, integrating from r = oo to r = with j3' e (k, oo) 



Pe(k, oo) = 0, to get the exact phase shifts. To get the Born approximation to /3e, we 
solve the equation iteratively, writing Pi = gfln + g 2 /3i2 + .-■■, where (3a satisfies 

- %0' tx - 2ik Pe (kr)/3' ei + ~ X (r) = (5.26) 



and Pi2 satisfies 



- i(3'; 2 - 2ik Pe {kr)(3' e2 + 2{f3' ll f = 0. (5.27) 



We can solve efficiently for fyi and (3^2 simultaneously by combining these two equa- 
tions into a coupled differential equation for the vector (fln, fin)- This method is 
much faster than calculating the Born terms directly as iterated integrals in r and 
will generalize easily to a theory requiring higher-order counterterms. 

Having found the phase shifts, we then use Levinson's theorem to count bound 
states. We then find the energies of these bound states by using a shooting method to 
solve the corresponding eigenvalue equation. We use the effective range approximation 
to calculate the phase shift and bound state energy near the threshold for forming 
an s-wave bound state. 



5.3 Results 

For the model at hand, we calculated the energy S[h] for a two-parameter (d and w) 
family of Gaussian backgrounds 

h(r) = -dve~ r2v2/2w2 . (5.28) 



In Fig. |5-2| , we show results which are representative of our findings in general. We 
plot the energy of configurations with fixed d — 1 as a function of w, for g = 1, 2, 4, 8 
(from top to bottom). We note that to this order, for g = 8 the vacuum is unstable 
to the formation of large = regions. 
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Figure 5-2: £[h) in units of v, for d — 1 and g = 1, 2, 4, 8, as a function of w. 

To explore whether the charged scalar forms a non-topological soliton in a given <p 
background, we add to £[h] the energy of a "valence" \1/ particle in the lowest bound 
state. We then compare this total energy to M, the energy of the \1/ particle in a flat 
background, to see if the soliton is favored. This is the scalar model analog of t-quark 
bag formation. 

For fixed g and m, we varied h looking for bound states with energy E such that 
E + £[h] < M. However, for those values of g and m where we did find such solutions, 
we always found that by increasing w, we could make £ [h] < 0, so that the vacuum is 



unstable, as we pointed out above in the case of g = 8 in Fig. [5-2[ Thus we find that 
if we stay in the g, m parameter region where the vacuum is stable, the minimum is 
at h = 0, so there are no nontopological solitons. 

Although we did not find a non-trivial solution at one-loop order in this simple 
model, our calculation demonstrates the practicality of our method in three dimen- 
sions. We can effectively characterize and search the space of field configurations, 
h(r), while holding the renormalized parameters of the theory fixed. The same meth- 
ods can be used to study solitons in theories with richer structure in three dimensions. 
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5.4 Derivative expansion 



Our results are exact to one-loop order. The derivative expansion, which is often 
applied to problems of this sort, should be accurate for slowly varying h(r). We 
found it useful to compare our results with the derivative expansion for two reasons: 
first, we can determine the range of validity (in d and w) of the derivative expansion; 
and second, where the derivative expansion is expected to be valid, it provides a check 
on the accuracy of our numerical work and C++ programming. Where expected, the 
two calculations did agree to the precision we specified (1 %). 

In our model, the first two terms in the derivative expansion of the one-loop 
effective Lagrangian can be calculated to be 



Ci = C ct + az + (3z + 



, , <fr ] f,. ,. 3 



32tt 2 



1 + z) 2 ln(l + z) - z - -z 2 



-1 l —( d z f 

48^ l + z l 11 ] ' 

(5.29) 



where z = gx/M 2 = (h 2 + 2hv)/v 2 , a and (3 are cutoff-dependent constants, and 
C ct is the same counterterm Lagrangian as we used in Sec. 2. For 4 scalar field 
theory a similar result was first derived in ||. The last term above is proportional 
to (dh) 2 , and is completely cancelled by a finite counterterm that implements the 
renormalization prescription of Sec. 2. In this prescription, counterterms also cancel 
the az and j3z 2 terms above. Thus the 0(p 2 ) derivative expansion for the effective 
energy, to be compared with the phase shift expression for S[h], is 



£de[Ji] 



d x 



9 2 v i 

327T 2 



m 



+ —,(h 2 + 2hvY 



[1 + z) 2 ln(l + z) 



z z 

2 



(5.30) 



The results of the comparison with the phase shift method can be seen in Fig. 



for d = 0.25, and g — 4. A similar pattern holds in general for other values of both 
d and g. As the width becomes larger, the two results merge. This is as expected, 
since it is for large widths, and thus small gradients, that we expect the derivative 
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Figure 5-3: \{£ — £de)/(£ — £ c i)| for d = 0.25, g = 4, as a function of w. 

expansion to yield accurate results. As the width tends to zero, both results tend 
to zero, and the fact that the plot tends to 1 simply indicates that the derivative 
expansion result goes to zero faster than the phase shift result. 
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Chapter 6 



Fermions in three dimensions 



It is now an easy matter to extend these three dimensional results to fermions. We 
continue to consider a classical background 0, now coupled to a Dirac fermion \1>. Our 
Lagrangian becomes 

c = \{d,<pf- ±(4? - v y + - 9<i»N> 

+a{d ll (b) 2 -b{(b 2 -v 2 )-c{(b 2 -v 2 ) 2 . (6.1) 

The small oscillations are now given by the Dirac equation, 

7° (-i^di + m + gh)^ k = ujip k (6.2) 

where h — (f> — v and m = gv. We will continue to work with a symmetric background 
h, so we can decompose our small oscillations into channels according to angular 
momentum. 

Letting a 1 = 7°7* and (3 = 7 , we will work in the basis 

/ (7fe \ / 1 , 

a k = and (3 = (6-3) 

\cj k Vo -1 
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where the a k are the standard Pauli matrices. We define the operator 



Q = (3{a- L) (6.4) 

where L = rxp. This operator commutes with the Hamiltonian, and is related to the 
total angular momentum operator J = L + ~<7 by Q 2 = J 2 + ~. We can thus restrict 
to states with a particular eigenvalue q of the operator Q, with q — ±1, ±2, ±3 
We also introduce the radial momentum operator 

Pr = -»•(-£■ + -) (6.5) 

V ar r / 



which allows us to rewrite eq. (|6.2j ) as 



(a r Pr + % -&rfiq + f3{m + gh(r)) \ ip k = u?p k (6.6) 



where a r — r ■ a. The solutions to this equation will be of the form 



'v - .... :; M w s 



r V ^(r)V A / 



for g > and 



for q < 0, where 3^f is the two-component spinor spherical harmonic corresponding 
to a state with total angular momentum j, orbital angular momentum £, and a z 
component of angular momentum M. These two solutions are eigenstates of parity 
with eigenvalues ±(— l) q respectively. The functions F and G obey the coupled radial 
equations 

((u-m)-gh(r))F+^ + ^-G = 

ar r 

((cu + m) +#a(r))G- — + = 0. (6.9) 

dr r 
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If we set h(r) = 0, F and G obey the second-order equations 

\ dr 2 r 2 J 

where u = \Jk 2 + m 2 . Thus the free solutions for F and G will be spherical harmonics 
ji(kr), with £ equal to the orbital angular momentum of the spherical harmonic that 
the radial function multiplies in eq. ( |6.7D or eq. (|6.8| ). 

In each channel, it is convenient to work with the full second-order equation that 
we obtain by squaring the Dirac equation. Any solution to the Dirac equation must 
solve this equation as well, and if we choose our boundary conditions so that the Dirac 
equation is solved, a solution to the second-order equation will automatically solve the 
Dirac equation. To facilitate our numerical computation, we will focus on the upper 
component, since its equation will involve the quantity u + m rather than uj — m, and 
the former is easier to deal with numerically at small k. (The small k limit is also 
the nonrelativistic limit, and in this limit the upper component is large compared 
to the lower component, which makes it easier to deal with numerically; if we were 
investigating negative energy solutions it would be easier to use the lower component.) 
We will therefore label our channels by parity (as above) and by I = 0, 1, 2, . . ., the 
orbital angular momentum of the upper component. The total spin j is then equal 
to £ + \ for positive parity and t — ~ for negative parity. (Only positive parity exists 
for t = 0.) 

Plugging F(r) = e^^rhg (At) into the full second-order equation for F gives 



, s +/ , 1 / 9 h'(r)(k Pi (kr) - ^±± + ipfY 



(6.11) 



in the positive parity channel, and 



f , K9 1 / , h(r)'(kpXkr) + l +i87') 



(6.12) 
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in the negative parity channel. As in the bosonic case, we have defined 



X = h 2 + 2hv. (6.13) 

Again, the phase shift is given by 

8f(k) = -2 Re f3f(k, 0). (6.14) 

We also have stronger divergences than in the bosonic case, so we must subtract the 
first four terms in the Born approximation and add them back in as diagrams. As 
a result, the numerical computation is somewhat more complicated; nonetheless, the 
basic method is the same: we iterate the differential equations for the phase shifts to 
obtain the Born approximations, which again cancel the divergences at large k and 
I (though for small k and £, the Born approximation may deviate widely from the 
exact phase shift). As in the boson model, the Born approximations to the phase 
shift go to zero at k — 0, while the exact phase shifts go to n times the number of 
bound states in that channel. And, as we would expect, there is an overall Fermi 
minus sign in front of the entire result. 

This computation adds the full Casimir energy to the computation of ||22|| ; as indi- 



cated above, to consistently compute the effect of the shift in the valence quark level 
requires that we include the Casimir energy as well. However, in this model, although 
the quantitative results are modified, the qualitative conclusions for moderate values 
of g are unchanged. 

An interesting generalization of this model is to add isospin, making the back- 
ground field into a complex 2-component vector. We can then construct "hedgehog" 
configurations that have nontrivial topology. These configurations are no longer CP 
invariant. In particular, they can have spectral asymmetry, which is reflected by a 
fermion level crossing during the process of adiabatically constructing the configura- 
tion from the trivial vacuum. If a level crossing does occur, the configuration carries 
fermion number and we do not need to fill an energy level at all. Even if a level 
crossing does not occur, the level that is heading toward zero is often strongly bound 
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and thus costs only a small amount of energy. 

To consider configurations that violate CP, we must address the question of 
whether to sum the quantity — ^\u\ over all modes or the quantity u the negative 
modes. (If C or CP is not broken, the two methods give the same answer.) The for- 
mer shows more consistency with the boson computation, while the latter is what one 
would expect from a Dirac sea picture, since the negative modes are the ones that are 
filled. The first method is correct. One can prove this formally |2T|, but here we will 
show it via a simple example, ordinary QED. Consider QED with a massive fermion, 
and imagine turning on a very weak, slowly varying electric field. To leading order 
in the coupling, this field will shift both positive and negative energy levels in the 
same direction, with corresponding particle and antiparticle states being shifted by 
the same amount. If we were to sum only the shifts in the negative energy levels, we 
would find a divergent result at the leading order in the Born approximation. How- 
ever, there is no available counterterm to cancel this infinity. If instead we calculate 
the shift summing all of the levels, to first order we simply get zero, since the shifts in 
the particle states cancel against the shift in the antiparticle states. This result agrees 
with the loop calculation, where we find that the leading contribution is identically 
zero, with no contribution from the counterterm. (Note that this result is different 
from what we find to leading order for scalar field. In that case, the particle and 
antiparticle levels move toward zero, and both sums give the same divergent result. 
However, in this case there is a counterterm available to cancel out this divergence.) 

A hedgehog configuration breaks rotational invariance, but we can choose it to 
stay within the spherical ansatz, so that it is invariant under grand spin (simultaneous 
rotations in physical space and isospin space). We decompose the spectrum into a 
sum over channels labeled by their grand spin and their parity. Each channel now 
has two coupled degrees of freedom instead of one, but still can be described by a 
second-order (matrix) differential equation. Summing the two eigenphases gives us 
the total contribution to the density of states from that channel. This matrix equa- 
tion must also be iterated through fourth order to obtain the Born approximations 
corresponding to the divergent diagrams. 
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Calculation in this generalized model in principle requires no further modifications 
to the framework presented above. However, there are practical difficulties associated 
with the increased computational complexity. Even in the model without isospin, the 
third- and fourth- order diagrams involve multidimensional integrals over external 
momenta and Feynman parameters. One possible simplification is to use a "toy" 
boson model with the same divergence structure as the third- and fourth-order dia- 
grams. If chosen properly, this model will exhibit the same local at second order as 
the fermions have at third and fourth order, so subtracting its second-order phase 
shifts and adding back in its second-order diagram (in the fermionic renormalization 
scheme, of course) can regulate the theory in a computationally simpler way. Such a 
manipulation is possible only because the divergences are simple, local functions of 
the background field; the finite parts of the bosonic diagram will of course be very dif- 
ferent from the finite parts of the fermionic diagrams. However, because we add back 
in everything we have subtracted, using counterterms fixed in the fermion theory, we 
have not changed the final answer. Work is underway on this model ||23|| . 
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Chapter 7 



Conclusions 



We have seen that the continuum density of states formalism has allowed us to apply 
the powerful calculational tools of quantum mechanics — phase shifts, Born approx- 
imations and Levinson's theorem — to nontrivial problems in quantum field theory. 
These formal tools provide the robust framework we need to do unambiguous cal- 
culations. Without them, we could have easily missed the subtleties of bosonic and 
fermionic spectra that we have seen above, and such a mistake would lead to a dras- 
tic change in the end result. Such errors are intolerable because of their theoretical 
consequences, as in the case of the saturation of the BPS bound in 1+1 dimensional 
super symmetric theories, and because of their phenomenological consequences, as in 
the case of heavy quark "bags" in the Standard Model. 

One can envision many other applications of this work. As mentioned above, 
the most immediate is to include fermions coupled to scalars with isospin in the 
spherical ansatz. One could also imagine adding gauge fields to such a model, to 
bring it still closer to the real Standard Model. One could then apply these methods 
to other gauge theory solitons, such as 't Hooft-Polyakov monopoles or flux tubes in 
superconductors. It would also be interesting to see if one could use these techniques 
to compute quantum corrections to the action of instantons. 
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Appendix A: Properties of 
reflectionless potentials in one 
dimension 



In this section we review the properties of solutions of 

/ d 2 ^ 



dx 2 



+ Ve(x) r](x) = k r](x) 



with 

V e (x) = -£{£ + l)m 2 sech 2 mx (2) 



for an integer £. For more details, see [19|. The results used in the soliton problems 



above will then be obtained by rescaling x — > |. The basic technique will be to define 
raising and lowering operators 

ai = i ( ^— + imtanhmx ] 
\dx J 

a\ = fmtanhmij (3) 

so that we can rewrite eq. ([[]) as 

(a\ae — £ 2 m 2 ^j r](x) = k 2 r](x) = {az + ia\ +1 — (£ + l) 2 m 2 ^j rj(x). (4) 

We can then connect the solutions for a given £ with those for £ — 1 by observing that 
eq. (f|) implies that if rj(x) is an eigenstate in the potential 14- i, then r](x) = a\fj(x) 
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is an eigenstate in the potential Vi with the same eigenvalue k 2 . Thus we have 



u e i](x) = a e r)(x) 

wevi x ) = a eV( x ) (5) 

with cut = \/k 2 + l 2 m 2 . Thus the spectra of Vg and V^_i are identical, except that 
the spectrum of Vg might contain additional states annihilated by ae. Indeed, there 
is exactly one such state, at k 2 = —£ 2 m 2 (which becomes the zero mode of the 
corresponding soliton). 

By iterating this process, we arrive at I = 0, which is a free particle. Since we 
know how to solve that problem completely, we can go backwards and solve for any Vg 
simply by applying raising operators and at each step solving for the one additional 
zero mode. 

For example, for i = 1 we start with the free scattering states e tkx and apply a\, 
giving 



rilix) = — ( mtanhmx | e tkx = (k + imta,nhmx)e lkx 

kK ' uji \dx ) uji ' 



(6) 



and a new bound state at k 2 = — 1 that satisfies 

i ^- — h mtanhmij Vo( x ) = ~^ Vo(k) = — ^=sechmx. 

We can then proceed to £ = 2, obtaining scattering states 

Vk( x ) = ~ ( j 2mtanhmx\ rjl(x) 

LO2 \dx J 

-((k + 2imtanhmx)(k + imtanhmx) + m 2 sech 2 mx)e ikx 



(7) 



(m 2 + k 2 + 3imkta,nhmx — 3m 2 tanh 2 mx)e lkx (8) 
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a bound state at k 2 = — 1 given by 

rjf(x) = —j= ^- 2mtanhm:rj i]o( x ) = —i]J^sechmx tanhmx (9) 

and a new bound state with k 2 = —4 solving 

i + 2mtanhmx^ i]q(x) = — > ^q(^) = y^sech 2 ma;. (10) 

We can extend to any integer £ by repeating this process. We find that the 
scattering from these potentials is reflectionless, with phase shift 

S e (k) = 2J2 tan- 1 ( 3 ^f) (11) 



j'=i 

and bound states at 



k 2 = -(mj£) 2 (12) 

for j — 0, . . . , £. The state at /c 2 = is a "half-bound" state, right at threshold, which 
goes to a constant as x — > ±oo. This state is given by acting with our raising operators 
on the threshold state in the free case, which is just a constant wavefunction. Thus 
for i = 1 we find 

^thresh = — mtanhmxj 1 = — itanhmx (13) 

and for £ = 2 we obtain 

^thresh (s) = ^ ^ - 2m tanh mx^J 77^^ (x) 

= ^(1 -3tanh 2 mx). (14) 

We see that in both cases, these states are the k — > limit of the scattering states, 
and have the property that they go to a constant (rather than a straight line of some 
nonzero slope) as x — > ±00. 
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Since Vi is symmetric, we can separate the spectrum of wavefunctions into sym- 
metric and antisymmetric channels. The symmetry of the bound states will alternate, 
with the lowest energy bound state being symmetric. The phase shift can also be de- 
composed into contributions from the two channels, with 

6i(k) = Sf(k) + Sf(k). (15) 



That the scattering is reflectionless is equivalent to 

5*(k)=5t(k). (16) 
Levinson's theorem for the two channels gives |7j 

5, A (0) = nnf 

<5?(0) = 7r(nf-i) (17) 

where nf and nf are the numbers of antisymmetric and symmetric bound states 
in partial wave £, with threshold states counted as one half. Thus we see that the 
threshold states are essential in reconciling eq. ( |T6| ) with eq. flT7p. 
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Appendix B: Phase shifts in 
arbitrary dimensions 



In this section we demonstrate explicitly the link between the leading Born approxi- 
mation to the phase shifts and tadpole diagrams by computing both in dimensional 
regularization. We consider only the Bose case for simplicity. 

In n spatial dimensions, we consider the bosonic small oscillations in partial wave 
£ with wave number k in the presence of a n— dimensionally spherically symmetric 
bosonic potential V(x). As in three dimensions, we can separate variables and obtain 
a differential for the wavefunction i](r) in partial waves t. Also as in three dimensions, 
it is convenient to use a rescaled wavefunction u = r^~ lr q, giving a generalized Bessel's 
equation 



This is the same equation as in three dimensions except that the eigenvalue of the 



The degeneracy factor for this partial wave (i.e. the dimension of this representa- 
tion of the rotation group, which gives the generalization of 2£+l in three dimensions) 
is the dimension of the space of symmetric tensors with t indices that each run from 
1 to n with all traces (contractions) removed. Working out the combinatorics gives 




u + Vu = k 2 u. 



(18) 



Casimir operator L 2 = ^M a pM alS has been generalized from 1(1 + 1) to £(£ + n — 2). 



N? = 



(n + £- 1)! 
£\(n-l)\ 



(£-2)\(n- 1)! 



(n + £-3)\ 



(19) 
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which we will analytically continue to arbitrary n using the gamma function, which 
satisfies 

xT (x) = T (x + 1) and T (n + 1) = n! (20) 

For a fuller derivation of eq. ([T8|) and eq. ([TP]) see for example Appendix B of [El] . 

As we have seen above, the tadpole graph requires the external momentum to be 
equal to zero. Thus we should expect that both the leading Born approximation and 
the tadpole graph will depend only on the spatial average of the potential (its p = 
component), 

2n 
7T 2" pOO 

{V) = - T - K V(r)r^dr. (21) 



T ( 2 ) Jo 



The leading Born approximation to the phase shift is 



CO 



5f\k) = ~ I J %+ t-i{krfV{r)rdr (22) 
2 jo 2 

where J v (z) is an ordinary Bessel function, related to the spherical Bessel function 
by jp( z ) — \fi^Jp+±( z )- The leading correction to the energy is 



em =EN ,r^. m) ^_m (23) 

Jo 2n dk 



where we have used Levinson's theorem as in eq. ( |2.10|) . This manipulation was not 



necessary in our calculations in three dimensions since we did not have any singular- 
ities at k — 0, but it still would have been a perfectly valid thing to do; all the Born 
approximations vanish at both k = and k = oo. However, to calculate in arbitrary 
dimensions we will need to perform the I sum first. Thus, to avoid the same infrared 
singularities we saw in one dimension, we again have used Levinson's theorem to add 
a total derivative to the unregulated integral. As a result, our analysis of arbitrary 
dimensions will provide further evidence that the correct form of the unregulated sum 
is indeed that of eq. ( |2.10|) . 
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Using the Bessel function identity 

^ (2g + 2l)(2g + l-l)! (2g)! /z\ 2 * 

fro ii ) "WUJ (24) 

we can explicitly do the sum over £, giving 

£ {1) = r ["(to -m)(n-2)k n - 3 dk. (25) 

(4vr)fr(f)io 

The fc integral can be calculated near n—\ and then analytically continued, giving 

f(— >^*~ S r (^) r (^) «> 

and so we find 

= :4r fi^) ™- (27) 



(47r) 2 V V 2 

in agreement with what we obtain using standard dimensional regularization of the 
tadpole diagram in n + 1 space-time dimensions. 



70 



Bibliography 



[1] E. D'Hoker and E. Farhi, Nucl. Phys. B248 (1984) 59, E. D'Hoker and E. Farhi, 
Nucl. Phys. B248 (1984) 77 

[2] J. Bagger and S. Naculich, Phys. Rev. Lett. 67 (1991) 2252, Phys. Rev. D45 
(1992) 1395. 

[3] J. Schwinger, Phys. Rev. 94 (1954) 1362. See also W. Greiner, B. Muller and 
J. Rafelski, Quantum Electrodynamics of Strong Fields (Springer, Berlin, 1985). 

[4] E. Farhi, N. Graham, P. Haagensen and R. L. Jaffe, Phys. Lett. B427 (1998) 334- 
342, N. Graham and R. L. Jaffe, Phys. Lett. B435 (1998) 145-151, N. Graham 
and R. L. Jaffe, Nucl. Phys. B544 (1999) 432-447, N. Graham and R. L. Jaffe, 
hep-th/9901023, to appear in Nuclear Physics B. 

[5] S. Coleman, Aspects of Symmetry (Cambridge University Press, Cambridge, 
1985); (North-Holland, Amsterdam, 1982). 

[6] L. Schiff, Quantum Mechanics (McGraw-Hill, NY, 1968). 

[7] G. Barton, J. Phys. A: Math. Gen. 18 (1985) 479. 

[8] A. Rebhan and P. van Nieuwenhuizen, Nucl. Phys. B508 (1997) 449. 

[9] H. Nastase, M. Stephanov, P. van Nieuwenhuizen and A. Rebhan, 
Nucl. Phys. B542 (1999) 471. 

[10] M. Bordag, J. Phys. A: Math. Gen. 28 (1995) 755. 



71 



[11] R. Dashen, B. Hasslacher, and A. Neveu, Phys. Rev. D10 (1974) 4114, 4130. 

[12] A. d'Adda and P. di Vecchia, Phys. Lett. 73B (1978) 162; A. d'Adda, 
R. Horsley and P. di Vecchia, Phys. Lett. 76B (1978) 298; R. Horsley, 
Nucl. Phys. B151 (1979) 399; S. Rouhani, Nucl. Phys. B182 (1981) 462.; 
H. Yamagishi, Phys. Lett. 147B (1984) 425; A. K. Chatterjee and P. Majum- 
dar, Phys. Rev. D30 (1984) 844; Phys. Lett. 159B (1985) 37; A. Uchiyama, 
Nucl. Phys. B244 (1984) 57; A. Uchiyama, Progr. Theor. Phys. 75 (1986) 1214; 
A. Uchiyama, Nucl. Phys. B278 (1986) 121; A. Rebhan and P. van Nieuwen- 
huizen, Nucl. Phys. B508 (1997) 449. 

[13] R. K. Kaul and R. Rajaraman, Phys. Lett. 131B (1983) 357; 

[14] J. F. Schonfeld, Nucl. Phys. B161 (1979) 125. 

[15] C. Imbimbo and S. Mukhi, Nucl. Phys. B247 (1984) 471. 

[16] D. Olive and E. Witten, Phys. Lett 78B (1978) 97. 

[17] K. Schoutens, Nucl. Phys. B344 (1990) 665; C. Ahn, D. Bernard and A. LeClair, 
Nucl. Phys. B246 (1990) 409; C. Ahn, Nucl. Phys. B354 (1991) 57. 

[18] R. Jackiw and C. Rebbi, Phys. Rev. D13 (1976) 3398. 

[19] P. Morse and H. Feshbach, Methods of Mathematical Physics (McGraw-Hill, 
New York, 1953). 

[20] M. Shifman, A. Vainshtein and M. Voloshin, Phys. Rev. D59 (1999) 045016. 

[21] R. Alkofer, H. Reinhardt, and H. Weigel, Phys. Rept. 265 (1996) 139. 

[22] S. Dimopoulos, B. Lynn, S. Selipsky and N. Tetradis, Phys. Lett. 253B (1991) 
237; 

[23] S. Bashinsky, E. Farhi, N. Graham, R. L. Jaffe and H. Weigel, to appear. 

[24] S. Lee, S. Minwalla, M. Rangamani, and N. Seiberg, Adv. Theor. Math. Phys. 2 
(1998) 697. 



72 



